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Abstract. In this paper we present a fully selfconsistent 
SSC model with particle acceleration due to shock and 
stochastic acceleration (Fermi-I and Fermi-II-Processes re- 
spectively) to model the quiescent spectral energy distribu- 
tion (SED) observed from PKS 2155. The simultaneous 
August/September 2008 multiwavelength data of H.E.S.S., 
Fermi, RXTE/SWIFT and ATOM give new constraints to 
the high-energy peak in the SED concerning its curvature. 
We find that, in our model, a monoenergetic injection of 
electrons at 70 = 910 into the model region, which are ac- 
celerated by Fermi-I- and Fermi-II-processes while suffering 
synchrotron and inverse Compton losses, finally leads to the 
observed SED of PKS 2155-30.4 shown in lH. E. S. S."andl 
|Fermi-LAT collaborations! ( |2009 V In contrast to other SSC 
models our parameters arise from the jet's microphysics and 
the spectrum is evolving selfconsistently from diffusion and 
acceleration. The 70-factor can be interpreted as two coun- 
terstreaming plasmas due to the motion of the blob at a bulk 
factor of r = 58 and opposed moving upstream electrons at 
moderate Lorentz factors with an average of 7„ « 8. 



1 INTRODUCTION 

Among the class of active galactic nuclei (AGN), blazars are 
showing a spectral energy distribution (SED) that is strongly 
dominated by nonthermal emission across a wide range of 
wavelengths, from radio waves to gamma rays, and rapid, 
large-amplitude variability. Presumably, these characteristics 
are due to a highly relativistic jet which covers a small angle 
to the line-of-sight, emitting the observable Doppler-boosted 
synchrotron and inverse Compton radiation. 
In SSC models the characteristic double humped spectra 
of blazars are explained by electrons in the jet emitting 
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synchrotron radiation while being accelerated in a magnetic 
field, which gives the first peak in the SED. These high 
energy electrons upscatter the very same synchrotron 
photons to TeV energies due to the inverse Compton effect, 
resulting in a second peak in the SED. Another approach to 
the explanation of the double humped structure are proton 
initiated electromagnetic cascades (e.g., Mannheim 1993| l 
or external Compton models. 

The key issue is to understand which physical mechanisms 
are leading to such SEDs. In particular that means to explain 
the ultrarelativistic electron spectra within the jet, which are 
believed to be responsible for the gamma radiation. 
The high peaked BL Lac objects (HBLs) as a subclass 
of blazars show a peak in their SED in the X-ray regime, 
suggesting that an inverse Compton peak should occur at 
correspondingly high gamma-ray energies. In fact, a large 
fraction of the known nearby HBLs have already been 
discovered with Cerenkov telescopes, such as H.E.S.S., 
MAGIC, and VERITAS. Since 2008 the Fermi satellite 
measures at these high gamma-ray energies. The energy 
range of the Fermi data is slightly different from the H.E.S.S. 
and VERITAS-Telescopes which gives new constraints to 
the SEDs. 

The first Fermi data pubUshed is from PKS 2155-30.4, 
a HBL at redshift z = 0.117 (luminosity distance: 
4 = 1.67 • lO^"^ cm) ( |H. E. S. S. and Fermi-LAT coUabora- 
[tioiisl [2009) 1. 

We present a selfconsistent SSC model that is not only 
able to model the SED of PKS 2155-30.4 shown in 
|H. E. S. S. and Fermi-LAT collaborationsl ( |2009] l but also to 
partly explain the "ad-hoc" injected particle spectra of many 
SSC models. Therefore we introduce and solve the kinetic 
equation describing the synchrotron-self-Compton emission 
numerically in two different zones within the jet (see section 
[2]). We use the exact Klein-Nishina cross section which 
is important at the relevant very high gamma-energies to 
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describe the inverse Compton radiation and energy losses of 
the electrons. The emphasis lies on the accurate treatment of 
the two possible particle acceleration mechanisms (Fermi-l- 
and Fermi-II) which are able to produce high energy 
electrons as well as on the selfconsistent treatment of the 
radiation processes. 

2 THE MODEL 
2.1 Model geometry 

We extend the well-established SSC model by Fermi-I and 
Fermi-II acceleration mechanisms to a selfconsistent SSC 
model with two zones in a nested setup. Both regions (the 
acceleration- and the radiation zone) forming the blob are as- 
sumed to be spherical and homogeneous containing isotrop- 
ically distributed non-thermal electrons and a randomly ori- 
ented magnetic field. The acceleration zone is assumed to 
be spatially significantly smaller than the surrounding radi- 
ation zone. Furthermore every electron leaving the acceler- 
ationzone enters the radiation zone. These assumptions are 
common place in SSC models (e.g., |Kirk et al.[[T998^ . 
To derive the kinetic equations describing the time evolution 
of He (7), Ne{'-f) (rie in the acceleration zone, in the ra- 
diation zone) as the differential electron densities we use the 
one dimensional diffusion approximation (eq. ([T]l) of the rel- 
ativistic Vlasov equation (e.g. ' Schlickeiser| |2002| l, which is 
applicable due to the assumptions made above. 
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where f{p, t) is a particle distribution function, with the par- 
ticle's momentum value p. F describes the contributing pro- 
cesses, such as synchrotron radiation or acceleration, in mo- 
mentum space. Catastrophic particle gains and losses are 
considered via S{p, t). 

Making use of the relativistic approximation E pc — -fmc 
and the relation n{p,t) — Anp"^ f{p,t) one can derive the 
kinetic equations governing the model. 

2.2 Kinetic equations 

2.2.1 Acceleration zone 

While the blob propagates through the jet, electrons are con- 
tinuously injected into the acceleration zone when consider- 
ing the blob's rest frame, leading to an injection function 

Qinj(7,0 := QoSh - loMt -to) ,to = (2) 

which we assume monoenergetic and time independent. 
These low- to mid-energy electrons are accelerated system- 
atically and stochastically due to Fermi-I and Fermi-II pro- 
cesses while suffering synchrotron and inverse Compton 
losses. Energy losses due to inverse Compton scattering are 



calculated using the full Klein-Nishina cross section, see eq. 
( fTS) . This leads to Pic (7) given in eq. (|8]l with the corre- 
sponding radiation field npu of the acceleration zone. Due 
to the non equilibrium of magnetic and radiative energy in 
the acceleration zone the energy losses via inverse Comp- 
ton scattering can become quite significant and must not be 
neglected, also the Thomson limit is not appropriate here. 
The synchrotron losses are calculated using eq. (|3j from 
|Grnzburg and Syrovatskii| ( |1969[ ) for isotropic particle dis- 
tributions 
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with the Thomson cross section uj. According to SchUck- 



eiser ( 1984 ) particle acceleration via parallel shockfronts and 
stochastic acceleration caused by scattering at Alfven waves 
leads to 
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for the function F. With the parallel spatial diffusion coeffi- 
cient i^i I , which is momentum independent for hard spheres 
and the characteristic speeds va for the Alfven mediated 
stochastic acceleration and vs for parallel shockfronts. Sub- 
stituting p — > 7 in eq. ([T} and eq. (|4]) according to the rela- 
tivistic approximation mentioned above, one will finally find 
eq. (|5]i; the kinetic equation of the acceleration zone. 

+ P^cil) ~ t.-h) • ne(7, t)] + 
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where the characteristic acceleration timescale t^^^c is given 
by 



(6) 
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Eq. (|6]l for t^^c is a direct consequence of the derivation of eq. 
Q out of eq. ([TJ using eqn. Q and The expression in 
eq. (|6]l includes the analytical timescale for non-relativistic 
shock acceleration. Ac cording to [ Bednarz and Ostrowski 
1996[ ) and especially to Ellison et al. (1990i the accelera 



tion timescale for parallel relativistic shock waves decreases 
approximately by a factor of 3. We did not take into account 
this behavior for it is unclear how the analytical expression 
looks like in that case. Secondly we are omitting the energy 
dependency of t^^c using hard spheres for the plasma insta- 
bilities anyway. This issue is irrelevant for the modelling 
(for we are setting numerical values for tacc) and the type of 
energy spectrum (a powerlaw) produced by Fermi-I acceler- 
ation is identical in the non-relativistic and relativistic case 
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( Sokolov et al. 2004 1. But it has to be kept in mind for the 
interpretation of the results. 

The parameter a w v1/v\ determines the ratio of shock to 
stochastic acceleration, tesc = V^acc/c is the characteristic 
timescale for electrons escaping from the acceleration region, 
where rj is an empirical factor set to 77 = 10 and i?acc the ra- 
dius of the acceleration sphere. All escaping electrons enter 
the radiation zone downstream the jet. The seperation in two 
zones can firstly explain the injected electron spectra and sec- 
ondly takes account of a much more confined shock region 
for Fermi-I acceleration will probably not occur in the whole 
blob region when considering physical sources. 
Our model can be compared with the model presented by 
Katarzynski et al.| ( |2006 1. The kinetic equation (eq. 3 in their 
paper) is almost similar to the kinetic equation in the acceler- 
ation zone eq. [5] One major difference to our model is their 
sole use of stochastic acceleration. In fact their model is the 
limit of our model for a — > 0. Additionally they limit them- 
selves to radiation in the acceleration zone, which is useful 
when not taking into account shock acceleration. Besides 
that there are number of minor differences regarding the ex- 
act treatment of inverse Compton losses and the derivation of 
escape rates. 

Due to the small spatial extent the acceleration zone does not 
contribute to the SED directly, i.e. riph(j^) is only calculated 
in order to determine the inverse Compton loss rate for the 
electrons in the acceleration zone. 

2.2.2 Radiation zone 

The electrons are not accelerated here. Thus the kinetic equa- 
tion takes the simple form 
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Electrons in the radiation zone suffer synchrotron (cx /3s7^) 
and inverse Compton losses (eq. ([S)), other energy losses are 
irrelevant in jetsystems of such low electron density (see e.g. 



Bottcher and Chiang 2002 1. 
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The integrals in eq. ([8]) are solved numerically using the 
full Klein-Nishina cross section for a single electron given 
in eq. ([TSj. The photon energies are rewritten in terms of 
the electrons rest mass, i.e. hiy = amc^ for the scattered 
photons and hv — aim(? for the target photons. The 
integration bounds of the outer integral in eq. ([HJ are a direct 
consequence of the kinematics. The non-trivial dependency 
of ^'ic(7) from iVph and thus of TVe from eq. (|8]l makes 
the numerical treatment of the kinetic equations inevitable 



lossrates for a single electron 

in the steady state photon distribution of PKS2155 




Fig. 1. Lossrates due to the inverse Compton effect for an electron 
of Lorentzfactor 7 in the model photon field of PKS 2155 (Thomp- 
son limit, red and used Klein-Nishina treatment, black) compared 
to the synchrotron losses, blue. Other losses like adiabatic cooling, 
pair production are irrelevant at typical SSC configurations. 



leading to a time resolved model. The loss rates for the 
electron distribution of PKS 2155 in the steady state are 
shown in Fig. 1, which indicates that in our case the inverse 
Compton losses would be slightly overestimated in the often 
used Thomson limit for all Lorentzfactors 7 because of the 
dependency on the photon field Pthom J dvyN^if^"^ and 
its special shape due to the modified injection at 70 = 910. 
This would not be the case for low energetic injected elec- 
trons and the resulting powerlaw-like photon distribution. 
For high Lorentzfactors 7 however the deviation of the 
Thompson limit for the inverse Compton scattering to the 
real Klein-Nishina treatment becomes more significant in 
each case. Again electrons escaping the blob are taken into 
account via fesc.rad — vRrad/c with the empmcal factor rj 
set to 77 = 10. Both electronic PDEs are connected via the 
catastrophic particle loss/gain-term. The factor (i?acc/^iad)'^ 
ensures particle conservation. 

To determine the time-dependent spectral energy distri- 
bution of blazars we solve the differential equation for 
the differential photon number density, obtained from the 
radiative transfer equation, including the corresponding 
terms with respect to the SSC model. 



= Rs- ca^iVph(z^, t) + Rc J . (9) 



dt 



^'ph.esc 



To describe the synchrotron photon production rate Rg in a 
converuent way we use the well known Melrose approxima- 
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tion (see e.g., [Pohll [2002 1 



R, ^ 1.8- 



Veil) 



(10) 



with the characteristic frequency eq. ( [TT| l for an electron with 
a Lorentz factor of 7. 
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In optically thick regimes the emitted synchrotron radiation 
is absorbed by the emitting electrons itself. This is described 
by the synchrotron self absorption coefficient, 



1 c , N,{jo,t) 
au = —^^PAlc) "2 • 
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(12) 



(with 7c = /(i^c)"^)- Here we made use of the monochro- 
matic approximation (e.g., |Felten and Morrison[ |1966| l for 
the synchrotron power: 

P.(.,7) = ^^^^ rd^'Kd^') (13) 



TOC I^cl Jo 

In SSC models the second hump in the SED of a blazar is due 
to inverse Compton scattered photons by the synchrotron ra- 
diation emitting electrons themselves. Here the full Klein- 
Nishina cross section from [Blumenthal and Gould] ( |1970| l 
is used to calculate the inverse Compton photon production 
rate. 
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To fully exploit the Klein-Nishina cross section, eq. ([14]), we 
used the approximate inverse Compton spectrum of a single 
electron scattering off a unit density photon field (e.g.,[Jonesj 
[T968i|Jauch and Rohrlich[[T976l l: 



rfiV(7,ai) 
dtda 



2Trrlc 
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2 (1 + 4ai79) 
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(15) 



with the electron's Lorentz radius tq — e^/(mc^), the 
scattering parameter q — a/(4Q!i7^(l — af-f)) and w 
1/(47^) < g < 1. Due to momentum and energy conserva- 
tion this equation is valid for ai < a < 4q!i7^/(1 + 4q;i7). 
The last catastrophic term in eq. (|9]l describes photons escap- 
ing from the emitting region, where 



t 



ph, esc 
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(16) 



is the approximate escape time, with _Riad the radius of the 
emitting blob. The escape time is chosen to be the light 
crossing time of the photons. 

The photon lossrate due to the pair production of electrons 
and positrons is not taken into account for two reasons. 
Firstly it is insignificant compared to the dominating 
synchrotron and inverse Compton processes. This is a 



consequence of the relatively low density (Bottcher and 



|Chiang[[2002[ ). Secondly it would violate the selfconsistency 
of our model for positrons are not treated, hence violating 
energy conservation. 

To compute the SEDs in our model we must shift the 
frame of reference from the blob to the observer For a 
sphere of radius R the observed flux at distance r is 



(17) 



With the Lorentz boosted intensity 7°'"* = (5^/^'°'' 
due to the bulk motion with a doppler factor S of the 
blob and the Lorentz transformed, red shifted frequency 
ly"'"' = ,5/(1 + z)v. Where 1^'°'' is calculated from the pho- 
ton unit density 



jblob 



hvc 
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for homogenous spheres. 
3 NUMERICS 

In our model we numerically solve the kinetic equations for- 
ward in time in order to obtain a model SED. The down- 
stream motion of the electrons induces the sequence of solv- 
ing the acceleration zone's equation before the kinetic equa- 
tion of the radiation zone in each time step. The simple Euler 
scheme was found adequate to do the time integration. 
In the acceleration zone we had to combine the Crank- 
Nicholson scheme ( |Press[ [2002[ ) with Godunov's method to 
provide both correct treatment of the characteristics and sta- 
bility for the derivation in 7. In the radiation zone the charac- 
teristic flows, due to the absence of acceleration, only in one 
direction making the Crank-Nicholson scheme sufficient. 
With our carefully tested code it is possible to calculate the 
dynamics of SEDs in a range of 20 orders of magnitude. The 
implemented code complies particle conservation in each 
zone alone and both together as well as the conservation of 
the total energy (i.e. of the electrons and the photons) over 
typical simulation times with a maximum error of 0(5%). 
For negligible stochastic acceleration, i.e. a 00, and with- 
out a radiation field, i.e. no inverse Compton losses, the 
steady state solution for the kinetic equation yields 



'^e,steady(7) = C* ^ ( A^ac 

T V7 
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with 7n]ax — (^acc/3s) and the constant C determined by 
the injection function Q. The implemented numeric model 
converges against this solution for sufficient simulation time. 
Additionally it was tested against the steady state analytical 
solution with Fermi-II processes given in |Schhckeiser| ( |7984l i 
with no significant deviations. Setting <esc,iad oo and ne- 
glecting inverse Compton scattering the spectral index of the 
powerlaw part of the electron distribution in the radiation 
zone is (analytically) reduced by one compared to the one 
in the acceleration zone, which also was confirmed by the 
implemented code. The inverse Compton scattering rate was 
confirmed against the approximate analytical results (low en- 
ergetic Thompson regime and extreme Klein-Nishina limit) 
before implementing. Concerning the photon distributions 
we validated the expected spectral indices in the steady state 
solution for the different frequency regimes, which together 
with the energy conservation between electrons and photons 
approves the integrity of the model. 

A detailed description of the used numeric techniques as well 
as the implemented model also in context with the variability 
of the sources will be given in a paper yet to be published. 



4 RESULTS 

The recent Fermi data give new constraints on the gamma- 
ray peak of the HBL PKS 2155-30.4 concerning its curva- 
ture. This is leading to a deep dip between the optical/X-ray 
and the gamma-ray peak. We are able to model the SED of 
PKS 2155-30.4 with our model by setting 



electron density 

acceleration zone 



70 = 910 



(19) 



for the monoenergetic injection into the acceleration zone. 
This is rather unusual but required to model the SED of PKS 
2155. Such moderate but not small Lorentz factors can be 
explained e.g. by two counterstreaming plasmas. If the up- 
stream electrons would be at rest, the bulk doppler factor of 
(5 = 116 would automatically lead to 70 = F « 58. As- 
suming speculatively that the upstream electrons moving in 
the opposite direction of the blob with a mean velocity of 
Vu hence a upstream Lorentz-Factor 7„ the 70 factor in the 
blob's rest frame must be calculated according to the rela- 
tivistic superposition: 



70 = 



\ 



1 - 



r2 + 72-i 



(20) 



Solving eq. pO| ) for our setup we find 7^ « 8 for the upsteam 
electrons which are streaming towards the blob. The numer- 
ically solved steady state electron density in the acceleration 
zone is shown in Fig. 2. We also show the time development 
for a "switched on" injection, i.e. 'T-e(7)lt<o = V-, and 
Q = Qo^il ~ 7o)i^(0' ™til the steady state is reached. 
In Fig. 2 it can clearly be seen tvhat accelerating electrons 
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injection function, 

time evolution 

steady state 




Fig. 2. Steady state electron distribution and its time development in 
the acceleration zone modelling the SED of PKS 2155-30.4 shown 
in Fig. 3 as arising from the injection function Qo- The corre- 
sponding intrinsic times are t = 1000 s (dashed black), t — 5000 s 
(dashed blue), t = 1 • 10* s (dashed red), t ^ 2 ■ 10^ s (dashed 
green). The steady state with its rising and falling powerlaw and the 
exponential cutoff at 7 ~ 10^ is reached at about t = 10^ s. 



using Fermi-I and Fermi-II processes leads to powerlaw elec- 
tron distributions with an exponential cut-off often used as 



the ad-hoc injection function in onezone-SSC models (Chi- 



ang and Bottcher 2002) l (right side of 70 in Fig. 2), thus ex- 



plaining them using the diffusion theory derived from plasma 
physics. By injecting electrons with eq. ^19\ and significant 
stochastic acceleration (i.e. a = 0(1)) we are also able to 
produce rising electron spectra before decreasing in a power- 
law and an exponential cut-off, like introduced in Bottcher] 
|and Chiang] ( |2002) (left side of 70 in Fig. 2). The Fermi-II 
processes are responsible for the rising power-law and expo- 
nential cut-off, whereas the ratio of iacc/^esc determines the 
spectral index of the power-law at 7 > 70. It can clearly 
be seen from Fig. 2 that the convergence against the steady 
state solution for the electron density begins relatively rapid 
while slowing down eventually. The simulation time when 
the steady state is reached corresponds to the escape time 
of the electrons in the acceleration zone. When concerning 
variability and time resolved lightcurves of blazars this is an 
advantage of the twozone model because the rising part in 
such lightcurves corresponds partially to the escape time of 
the acceleration zone (while the falling part is connected to 
the response time of the system, trad.esc)- 
An acceleration zone electron density as shown by the solid 
black curve in Fig. 2, leads to the desired broken power-law 
electron spectrum in the radiation zone which finally is able 
to model the SED of PKS 2155-30.4 (see Fig. 3 and Fig. 4). 
We used the parameters in Table[T]for the model SED in Fig. 
4 (black, solid line). The black dashed curve in Fig. 4 corre- 
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Table 1. Chosen parameters for the model SED shown in Fig. 4 
to fit the data ( |H. E. S. S. and Fermi-LAT collaborations! [2009l > of 
PKS 2155-30.4. 



]o{cm ^) B{G) 7?acc(cm) 7?rad(cm) tacc/W a r 



jection of electrons into the acceleration zone at time to = 
is shown in Fig. 3. It can clearly be seen that the final state 
of the model SSC correlates with the response time of the 
radiation zone tesc.rad and that the convergence again begins 
fast and slows down rapidly at higher simulation times. 



5.25-10'' 0.29 3.0-10" 6.3-10" 1.55 1 58 



time evolution (model) 

PKS2155 spectral energy distribution 




Fig. 3. Time evolution due to the switched on injection function 
until the steady state of the SED of PKS2155 is reached (see also 
Fig. 4). The intrinsic times are t = 1 ■ 10^ s (dashed black), t — 2 ■ 
s (dashed blue), t = 5 - lO"* s (dashed red), f = 1 - 10^ s (dashed 
green), the complete steady state is reached at about t — 2 ■ 10^ s, 
which correlates to the response time of the system due to tcsc.rad. 



sponds to a fit assuming a black body for the thermal contri- 
bution of the host galaxy thus the ATOM optical data is not 
to be taken into account for the SSC modelling. The curva- 
ture and deep dip in the model SED is a direct consequence 
of the rising part in the electron density of the acceleration 
zone. Thus it can be modeled by varying the ratio a of shock 
to stochastic acceleration. All the parameters in Table [T] are 
consistent with the limits given via other observations and 
statistics, e.g. determination of F using superluminal motion 
of Quasar jets. 

The recent Fermi, H.E.S.S. and ATOM data ( |H. E. S. S. and] 
|Fermi-LAT collaborations! !2009| l have been averaged over 
a period of 14 days and show a lowstate of the HBL PKS 
2155-30.4. This is confirmed by the | Aharonian et aL ( 2005[ l 
data of H.E.S.S. a few years ago which show the same flux 
level as the recent data. We used the EBL studies described 
in !Primack et al.! ( [2005| l to do the EBL deabsorption for the 
H.E.S.S. datapoints, a correction of the Fermi data is not nec- 
essary. 

The time development of the SED due to a switched on in- 



5 CONCLUSIONS 

Our model is able to explain the injection function of 
many onezone-SSC models as shock and stochastic ac- 
celeration of electrons upstream the jet entering the blob 
while continuously suffering synchrotron losses. By in- 
troducing Fermi-II acceleration we get rid of the sharp 
cut-off introduced in !Kardashevl ( |1962 |) or'Kirk et aL] ( |1998| l 
which probably does not occur in physical sources. Addi- 
tionally we are able to model relatively complex electron 
densities with increasing and decreasing parts through 
the stochastic acceleration of electrons, only by varying 
the monoenergetic injection to higher 70. In contrast to 
the ad-hoc injection of some onezone-SSC models such 
Lorentz factors have a physically reasonable, but highly 
speculative, explanation as upstream previously accelerated 
but already partially cooled electrons. These electrons are 
averagely moving in the opposite direction of the blob with 
a mean Lorentz factor of 7„ w 8 resulting, together with 
the motion of the blob, in 70 « 900 for the monoenergetic 
injection function used in the acceleration zone of our model. 

As recent data points out, these complex electron dis- 
tributions are necessary to model the new constraints 
concerning the gamma-ray peak of blazar's SEDs if one 
does not simply shift the synchrotron-peak to achieve the 
inverse Compton spectrum (e.g., !Kataoka et aL 2000) l. The 
curvature of the peak, and thus the deep dip between the 
two humps, is a direct consequence of the rising part in 
the responsible electron distribution within the blob. This 
constraint rules out many SSC models, which are not able to 
produce such electron spectra. 

With our model we are able to form the curvature of the 
gamma-ray peak and the dip by varying the influence of 
the Fermi-II processes. The shape and position of the 
synchrotron peak in the model SED is dominated by t^^.^. 
and i?acc, ^^rad as well as B. For the parameters concerning 
the acceleration arise from plasmaphysics considerations we 
gain insight into the jets microphysics while modelling ob- 
served SEDs. We have also shown that in such environments 
the Thomson approximation for the inverse Compton effect 
can not always be applied, especially when considering time 
resolution and hence non equilibria of the energy distribution 
in the blob. 

Here we only introduced steady state solutions of our 
model, but due to the spatially relatively small acceleration 
region, which is at least an oder of magnitude smaller than 
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the emitting region, this twozone-SSC model is able to 
selfconsistently model the rising part in the lightcurves of 
flaring blazars which are connected to the behavior in the 
acceleration zone, especially the energy transport from low 
to high energies. This, together with the consequences of 
the model geometry on the observable SEDs and lightcurves 
of blazars like in Sokolov et al. ( 2004| l, will be subject of a 
following paper. 
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Fig. 4. Lowstate of PKS 2155-30.4 with the simultaneous data of ATOM, SWIFT, RXTE, Fermi and H.E.S.S of the August/September 
2008 campaign from |H. E. S. S. and Fermi-LAT collaborations! ( |2009[ l (red triangles and circles). The 2003 H.E.S.S. data (blue circles) is 
also shown, proofing the lowstate of PKS2155-30.4. The VHE data have been deabsorbed using |Primack et aL] ( |2005^ . The dashed black 
curve shows a thermal fit for the contribution of the host-galaxy. Our model SSC fit, arising from the steady state electron distribution in the 
radiation zone is shown in the solid black curve, a moderate energy injection at 70 ~ 910 into the acceleration zone together with stochastic 
and systematic acceleration is needed to meet the curvature of the VHE peak given via the Fermi data. 



